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Abstract 

The detection of continuous gravitational-wave signals requires to ac- 
count for the motion of the detector with respect to the solar system 
barycenter in the data analysis. In order to search efficiently for such 
signals by means of the fast Fourier transform the data needs to be trans- 
formed from the topocentric time to the barycentric time by means of 
resampling. The resampled data form a non-stationary random process. 
In this communication we prove that this non-stationary random process 
is mathematically well defined, and show that generalizations of the fun- 
damental results for stationary processes, like Wiener-Khintchine theorem 
and Cramer representation, exist. 



1 Introduction 

Continuous gravitational-wave signals originating for example from spinning 
neutron stars are believed to be very weak. In order to detect them in the data 
of both bar and laser interferometric gravitational-wave detectors one needs 
to integrate the data for many days in order to achieve a signal-to-noise ratio 
sufficiently high to guarantee their detection. Over such time scales, however, 
the motion of the detector relative to the gravitational wave source introduces 
a Doppler modulation in the received signal. In order to effectively preserve 
coherent integration of the signal received by the detector, it is necessary in 
the data analysis to model very accurately its motion with respect to the solar 
system barycenter (SSB). 
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The simplest model of a monochromatic gravitational-wave signal emitted 
by such a source has the following form 

27r 

s = Acos[(f>o + 27r/ot + —foUo ■ rssB(i)] (1) 

In Eq. (|l|) (f>Q is the initial phase of the waveform, rggB is the vector joining 
the solar system barycenter (SSB) with the detector, np is the constant unit 
vector in the direction from the SSB to the neutron star. We assume that the 
gravitational wave form is monochromatic with frequency /o which we define 
as the instantaneous frequency evaluated at the SSB at i = 0. In general 
the frequency of the signal may change as a consequence of the spindown of the 
rotating neutron star. However, this does not introduce new qualitative features 
into the mathematical model we discuss here. 

The above signal cannot be detected by means of the fast Fourier transform 
algorithm (FFT) because the unknown frequency fo to be estimated is con- 
tained in the highly non-linear function of time depending on the position of 
the detector with respect to the SSB. Consequently the following technique is 
proposedo. Introduce a new time ts called barycentric time which is related 
to the topocentric time t by 

tB ^t + no-rssB{t)/c. (2) 

Then the signal given above takes the form 

SB = Acos{^a + 2T:fotB). (3) 

The signal ss is monochromatic and can be searched for by means of FFT. How- 
ever, as we shall show in the following section, the transformation given by Eq. 
H makes the noise in the data a non-stationary random process. The main re- 
sult of this communication, proven in the next section, will be to show that this 
non-sationary random process is mathematically well defined and generaliza- 
tions of the fundamental results for stationary processes like Wiener-Khintchine 
theorem and cramer representation exist. We also show that, if the spectrum 
of the unresampled noise is white, the resampled process remains stationary. 

Another case when the resampling of the data is applied are the so-called 
"accelarated searches" in pulsar data analysis S. These are searches for pul- 
sars in binary systems where the observation time is much shorter than the 
period of the orbit so that we can expand the phase of the signal in Taylor 
series keeping only the second order terms. This results in the following signal 

Sa^ Acos[(l)o+2TT{fot+Pot^)]. (4) 

This signal can be detected by matched filtering i.e. by multiplication of the 
data by exp{—2iTpt^) for a grid of prameters p followed by FFT. Nevertheless 
pulsar astronomers apply resampling also in this case. One defines a new time 

ta=t + at^, (5) 
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where a — Po/ fo- Then the signal becomes monochromatic: 



Sa = AcOs[(f>o + 271 fot a] 



(6) 



and can be searched for by means of FFT. Of course we do not know the 
frequency fo and parameter po, and therefore have to resample the data over 
an appropriate grid on the acceleration parameter a. 

2 Resampled random processes 

Let x{t) be a continuous parameter real-valued random process. The process 
x{t) is said to be completely stationary^ (sometimes referred to as strongly sta- 
tionary or strictly stationary) if for all n > 1, for any ti t2 i„ contained in the 
index set, and for any r such that ti + t, ^2 + ■■■,tn + t are also contained 
in the index set, the joint cumulative probability distribution function(cpdf) 
of x{ti), x{t2), x{tn) is the same as that of x(ti + t), x{t2 + t), x{tn + t). 
In other words, the probabilistic structure of a completely stationary process 
is invariant under a shift of time. The process x{t) is said to be second-order 
stationary«s (sometimes referred to as weakly stationary or covariance station- 
ary) if for all n > 1, for any tit2 t„contained in the index set, and for any 

T such that ti -\- T,t2 -\- ...,t„ + T are also contained in the index set, the 
joint moments of orders 1 and 2 of x{ti), x{t2), ...,x{tn) exist, are finite and are 
equal to corresponding joint moments of x{ti + T),x(t2 + t), ...,x(tn + t). A 
weakly stationary Gaussian random process is also completely stationary be- 
cause a Gaussian process is completely determined by its 1st and 2nd moment. 
For stationary random processes we have the following fundamental result. O 

Theorem 1 (The Wiener-Khintchine Theorem) A necessary and sufficient 
condition for p{t) to he the autocorrelation function of some stochastically con- 
tinuous (i.e. continuous in the mean square sense) stationary process, x{t), is 
that there exists a function, F{uj), having the properties of a distribution func- 
tion on (—00,00), (i.e. F(—oo) = 0, F(oo) = l,and Fiuj) non-decreasing), such 
that, for all t,p{t) may be expressed in the form, 



The necessary part of the above Theorem follows from a general theorem 
due to Bochner that any positive semi-definite function which is continuous 
everywhere must have the representation of the above form. In the case of 
a purely continuous spectrum we have dF{oj) = S{uj)dui, where S{uj) is the 
(normalized) spectral density function. Thus the Wiener-Khintchine theorem 
asserts that a well defined spectral density exists. 

Another important result is the existence of the spectral decomposition of 
the stationary random process itself. 




(7) 
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Theorem 2 (The Cramer representation) Letx{t),—oo < t < oo,6e a zero- 
mean stochastically continuous stationary process. Then there exists an orthog- 
onal process, Z{u)), such that for all t, x{t) may be written in the form, 

/■QC 
e"^*dZ{u), (8) 
-oo 

the integral being defined in the mean-square sense. The process Z{u)) has the 

following properties; 

(i) E[dZ{iv)] = 0, for all u 

(ii) E[\dZ{u)\'^] = dH{uj), for all u, 

where H{u}) is the (non-normalized) integrated spectrum of x{t), 
(Hi) for any two distinct frequencies, lu , uj' , [lj ^ lu'), 

cov[dZ{uj), dZ{(jj')] = E[dZ*{uj)dZ{u')] = 0. (9) 

The above result says that any stationary random process can be decomposed 
into a sum of sine and cosine functions with uncorrelated coefficients. 

Let x{t) be a random process and let tr = t-\-k(t\ 9) be a smooth one-to-one 
function both of the index t and the parameter set 6. Aas the random process 
x{t) taken at time tr i-e. y{tr) = x{t). In the following for simplicity we assume 
that the process x{t) is zero- mean. It immediately follows that the resampled 
process y{tr) is also zero mean. Suppose that the original process is stationary. 
Let us first convince ourselves that the resampled process is, in general, non- 
stationary. Let C{tr,tr) := E[y{tj.)y{tr)] be the autocovariance function of 
the resampled process. By definition of the resampled process we have that 
C{tj., tr) := E[x{t )x{t)] and by stationarity of x{t) we have C(i^, tr) = R{t —t). 
By implicit function theorem we have that there exists a smooth function t = 
tr + g{tr', 0). In order that the resampled process be stationary the function R 
must depend only on the difference {tr — tr) This is the case if and only if t 
is a linear function of tr i.e. t = tr + a{9)tr + b{9). Thus when the resampling 
transformation is non-linear the resulting resampled process is non-stationary. 

In the linear resampling case the Fourier transform Y{u)r) of the resampled 
process at frequency ujr is related to the Fourier transform X{(jj) of the original 
process at frequency w by the following formula: 

~, , expiw6~, . 

^K) = -f^x{uj), (10) 

where uj = -rr- 

l+a 

Let us consider the covariance function C{tr,tr), of the resampled process 
y{tr). It can be written as 

/oc 
<Pl,^{u;)<Pt^{u:)dH{w), (11) 
-oo ^ 

where we have introduced a set functions 

4>tr ('^) = exp[iLo{tr + g{tr))] (12) 
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We have the foUowing important result. 



Theorem 3 (General orthogonal expansions) Lety{tr) be a continuous pa- 
rameter zero mean process (not necessarily stationary) with covariance function 
C{tr,Sr) = E[y{tr)y{sr)]- If there exists a family of functions, {(f)t^.{^)}' de- 
fined on the real line, and indexed by suffix tr, and a measure, ^^.{lij), on the real 
line such that for each t, (j)t^{uj) is quadratically integrable with respect to the 
measure fi, i.e. 

/oc 
\4>t,Xuj)\^ dn{uu) < oo (13) 
-OO 

and for alltr, Sr^ C{tr^Sr) admits a representation of the form, 

/oc 
4>l,^{Lu)(pt^{uj)dH{u) (14) 
-OO ^ 

then the process y{tr) admits a representation of the form 

poo 

y{tr) - / cl>MdZ{co), (15) 



where Z(lo) is an orthogonal process with 

i;[|dZ(cj)p] = d/i(cj). (16) 

Conversely if yit^) admits a representation of the form with an orthogonal 
process satisfying , then C(t^, tr) admits a representation of the form . 



The formula (14) is called generalized Wiener- Khintchine relation and for- 
mula (|l5|) is called generalized Cramer representation of the random process. In 
our case the generalized Cramer representaion reads 

/•oo 

y{tr) ^ / eyip[iui{tr + g{tr))]dZ{uj), (17) 



This representation also clearly shows that the resampled process is in general 
non-stationary because the choice of basis function 

(ptri^) — exp{iu!tr) (18) 

is not in general possible. The generalized Cramer representation is also imme- 
diate from the Cramer representation for the original stationary process and its 
transformation to resampled time. It also follows that the measure ^ coincides 
with the integrated spectrum H{lj) of the original stationary process. However 
H{uj) cannot be interepreted as the spectrum of the resampled process. Indeed 
for the resampled process which is non-stationary the concept of spectrum is 
not mathematically well defined. 

The general orthogonal expansion theorem has already been used by M. B. 
Priestley [|| to develop the theory of so called evolutionary spectra. This theory 
describes a very important class of non-stationary processes often occurring in 
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practice for which the amplitude of the Fourier transform slowly changes with 
time. 

In the case of a continuous spectrum we have dH (lu) — S (uj) dio. Then one 
can write the properties (ii) and (iii) of the Cramer representation theorem as 



E 



[X{uj')*X{uj)]^ S{uj)S (^J -u'^ (19) 



where S is the Dirac function. In the continuous case it is instructive to calculate 
the correlation function for the Fourier frequency components of the resampled 
(non-stationary) process. Using Eq.(p^) the Fourier transform Y{lu) of the re- 
sampled process y(tr) can be written as 



Yiu) - / Qiui,iu)dZ{iOi), (20) 
where the kernel Q is given by 

Q{u!i,uj) = j (pt^ii^i) exp{—iujtr)dtr. (21) 



The correlation between two Fourier components of the resampled process takes 
the form 



oc /'OO 



E[Y{LOyY{Lu)] = I / Q{L0l,U;'yQ{LU2,L0)E[dZ{ui)*dZ{L02)] 

— oo J —oo 

Q{uji,Ljj'yQ{uJi,Ljj)S{uJi)duJi (22) 



Thus we see that for a resampled random process the Fourier components at 
different frequencies are correlated, This is another manifestation of the non- 
stationarity of the process. Let us now consider the example of white noise for 
which the spectral density is independent of frequency uj. It is straightforward 
to show that 

/>OC 

Q(wi, w')*Q(wi, w)dwi = 6{u' - u) (23) 



Thus for the case of white noise we have that 

E[Y{JyY{uj)]= S{u)6{lo' -u}j (24) 

and consequently in this case the noise remains stationary after resampling. It 
is also easy to see that the Fourier components at different frequencies will be 
uncorrelated if the spectral density is constant over the bandwidth of the kernel 
Q. It is possible that the last assumption will be fulfilled in the case of search 
of the data from gravitational-wave detectors from continuous sources^^ but is 
not guaranteed. 



6 



Acknowledgments 



We would like to thank Prof. Bernard F. Schutz for illuminating discussions 
concerning nonstationarity of the resampled noise. This work was supported in 
part by the KBN Grant No. 2 P03B 094 17. This research was performed at the 
Jet Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration (M.T.). 

References 

[1] p. R. Brady, T. Creighton, C. Cutler, and B. F. Schutz, Phys. Rev. D 57, 
2101 (1998). 

[2] H. M. Johnston and S. R. Kulkarni, Ap. J. 368, 504 (1991). 

[3] D. R. Lorimer, Radio pulsars - an observer's perspective in "The Neutron 
Star - Black Hole Connection" , NATO Advanced Study Institute, June 7-8 
1999, Elounda, Crete. 

[4] D. B. Percival and A. T. Walden, Spectral Analysis for Physical Applica- 
tions, Cambridge University Press, Cambridge 1993. 

[5] M. B. Priestley, Spectral Analysis and Time Series, Academic Press, Lon- 
don 1996. 

[6] See discussion in Section IV of Rcf. [Q . 



7 



